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Abstract. To explore the ionization conditions in highly-ionized absorbers at high red- 
shift we have studied in detail two intervening Ovi absorbers at z ~ 2 towards the quasar 
PKS 1448-232, based on high {R ^ 75, 000) and intermediate (R ^ 45, 000) resolution optical 
VLT/UVES spectra. We find that both absorption systems are composed of several narrow 
subcomponents with Civ/Ovi Doppler-parameters 6 < 10 kms"'^, typically. This implies 
that the gas temperatures are T < 10^ K and that the absorbers are photoionized by the UV 
background. The system at z = 2.1098 represents a simple, isolated Ovi absorber that has 
only two absorption components and that is relatively metal-rich {Z ~ 0.6 solar). loinization 
modeling implies that the system is photoionized with Ovi, Civ, and Hi coexisting in the 
same gas phase. The second system at z = 2.1660 represents a complicated, multi-component 
absorption system with eight Ovi components spanning almost 300 kms~^ in radial velocity. 
The photoionization modeling implies that the metallicity is non-uniform and relatively low 
(< 0.1 solar) and that the Ovi absorption must arise in a gas phase different from that traced 
by Civ, Cm, and Hi. Our detailed study of the two Ovi systems towards PKS 1448— 232 
shows that multi-phase, multi-component high-ion absorbers like the one at z = 2.1660 re- 
quire a detailed ionization modeling of the various subcomponents to obtain reliable results 
on the physical conditions and metal-abundances in the gas. 

Key words, galaxies: intergalactic medium - galaxies: quasars: absorption lines - cosmology: 
observations 



1. Introduction 

Highly ionized species like Ovi and CiV, ob- 
served in the spectra of distant quasars, are ex- 
cellent tracers of metal-enriched ionized gas in 
the filamentary intergalactic medium (IGM) and 
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in the circumgalactic environment of galaxies. 
Therefore, the analysis of intervening O VI and 
Civ absorbers towards low- and high-redshift 
QSOs is crucial for a better understanding of 
the physical nature, distribution, evolution, and 
baryon and metal content of the IGM in the con- 
text of galaxy evolution. Because of the high 
cosmic abundance of oxygen, the large oscilla- 
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tor strength of the O VI doublet (located in far- 
ultraviolet at AA1031.9, 1037.6 A), and the high 
ionization energies of the ionization states O"*"^ 
(113.9 eV) and 0+^ (138.1 eV), the Ovi ion 
is a particularly powerful tracer of the metal- 
enriched IGM and the gaseous environment of 
galaxies. Using QSO absorption spectroscopy, 
O VI absorption now is commonly detected in 
various different galactic and intergalactic envi- 
ronments in the redshift range z « — 3. 

In the local Universe, Ovi absorption in 
interstellar and intergalactic gas can be observed 
in the FUV spectra of stars and extragalactic 
background sources. For instance, Ovi absorp- 
tion is known to arise in the thick disk of the 
Milky Way (e.g. 



Savage et al.l 120031 ). in the 



extended, multi-phase gas halo s of the Milky 
Way and other galaxies (e.g . ISembach et al 
2003[ IWakker fc Savagd l2009t IProchaska et al 
20111 ). and in intervening Ovi absorption- line 
systems that trace meta l- enrich e d gas in the 



Tripp et al I I2OOOI: ISavage et all 
2OO2I: iRichter et al.l I200I ISembach et al.l \200^ 



IGM 



Tripp et al.l 120081: iThom fc Cheni " 12008c 
Danforth fc Shulll J2008[ ). for a review see 
Richter et al.l ( 2008f ). Over the last decade, 
intervening Ovi absorbers at low redshift were 
considered as major baryon reservoir in the IGM, 
possibly tracing shock-heated and collisionally 
ionized intergalactic gas tha t results from 
large - scale structure for mation ( Gen fc Ostrikeij 
I999I: iDave et al.l I2OOII) . This so-called warm- 
hot intergalactic medium (WHIM) has gas 
temperatures in the range 10^ < T < 10^ 
K and is believ ed to host 30 — 40% of the 
baryons at 2; = ( Gen fc Ostrikei 19991 ). Recent 
observational and theoretical studies indicate, 
however, that part of the O vi absorbers at 
low z may trace low-density, photoionized gas 
or conductive, turbulent, or shocked boundary 
layers between cold/warm 10'^ — IC* K) gas 
clouds and an ambient hot (~ 10^ — 10^ K) 
plasma ra ther than the shock-heated WHIM 
(|FoxI I2OIII see discussion in). Thus, a simple 
estimate of the ionization state of the gas 
in the absorbers from the observed Ovi/Hi 



ratios may lead to erroneous results because of 
the complex multi-phase character of the gas 
(jTepper-Garcia et al.ll201l[ ). 

For redshifts z > 2 Ovi absorption is de- 
tectable from the ground, where it can be ob- 
served in optical QSO absorption spectra at 
relatively high signal-to- noise (S/N). One very 
problematic aspect for the analysis of O vi ab- 
sorbers at high redshift is the often severe blend- 
ing of the Ovi absorption with the Lya for- 
est. As for low redshifts, the origin and na- 
ture of Ovi absorbers at high z is expected to 
be manifold. It has been shown by simulations 
e.g. Theuns et al. 2002; lOppenheimer fc Dave 



20081 ) that shock-heating by collapsing large- 



scale structures is not efficient at high redshift 
to provide a widespread warm-hot intergalac- 
tic phase in the early Universe. Instead, galac- 
tic winds probably contribute substantially to 
the population of photoionized and collision- 
ally ionized Ovi absorbers at high redshifts, 
enriching the surrounding circumgalactic and 
intergalactic gas with heavy el ements at rel- 
atively high gas tem peratures ([Fangano et al 



veiy 

l2007t Ikawata fc Ranch ,2007, ). In fact, many of 
the strong Ovi absorbers at high z exhibit 
complex absorption patterns that would be ex- 
pected for a cir cumgalactic multi-phase gas en - 
vironment (e.g. Bergeron fc Herbert-FortI 2005[) . 
Similar as for circumgalactic absorbers in the lo- 
cal Universe, a considerable fraction of the O vi 
absorbers at high z thus may arise in conductive, 
turbulent, or shocked boundary layers. 

Next to those Ovi absorbers that trace 
highly- ionized gas in the immediate environ- 
ment of galaxies, intergalactic O VI absorbers 
(i.e., absorbers that are not gravitationally 
bound to individual galaxies) may arise in 
regions that are sufficiently enriched with 
heavy elements. Previo us surveys of high- 
redshift O VI abso r bers ( Bergeron et al"] 2002t 
'Si mcoe et all l2002i [200I I2OO6I: ICarsweh et al ' 



12OO2I: iBergeron fc Herbert-FortI 120051 ) have 
shown that there are many narrow O VI 
absorbers with Doppler-parameters & < 10 
kms~^. Such narrow lines cannot arise from 
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coUisionally ionized gas but must be related to 
photoionized (possibly intergalactic) gas with 
temperatures T < 10^ K. Many of these narrow 
Ovi absorbers at low and high redshift display 
velocity-centroid offsets between O vi, C iv, and 
Hi, suggesting that these ions do not arise in 
the same gas phase. Unfortunately, this crucial 
aspect is only partially considered in previous 
O VI surveys. 

To explore the multi-phase character of high- 
ion absorbers and to improve our understanding 
of the ionization conditions in O vi systems it 
is important to investigate in detail the absorp- 
tion characteristics and ionization conditions in 
selected absorption-line systems. For this pur- 
pose, absorbers that can be observed at high S/N 
and for which the O vi absorption is not blended 
by Lya forest lines are particularly important. 
Because of the complexity of many high-ion ab- 
sorbers that often are composed of several veloc- 
ity subcomponents a spectral resolution of i? w 
45, 000 and higher is desired. Note that while 
the analysis of individual high-ion absorbers is a 
common strategy to explore t he nature of O VI 
absorbers at low redshif t (e.g. iTumlinson et al. 
20111 : ISavage et al. I l201lh . detailed studies of in- 
dividual O VI absorptioii syste ms at high redshift 



are rare fe.g.. lFox et al.ll2011l) 



In this paper we present VLT/UVES obser- 
vations at intermediate (i? « 45,000) and high 
(i? « 75, 000) spectral resolution of two particu- 
larly interesting O vi systems at z w 2 along the 
line of sight towards the quasar PKS 1448—232. 
This sightline was selected by us for a detailed 
study, as it contains two unsaturated Ovi sys- 
tems at Zabs = 2.1098 and Zabs — 2.1660, both 
displaying a well-defined subcomponent struc- 
ture with narrow O vi/C iv absorption compo- 
nents an d without major blend i ng with Lya for - 
est lines (|Bergeron et al.ll2002t IFox et al.il2008l ). 
These two absorption systems therefore repre- 
sent ideal targets to study in detail the physical 
conditions in photoionized, multi-phase high-ion 
absorbers at high redshift. 



2. Observations and absorption-line 
analysis 

2.1. VLT/UVES observations 

Our data set consists of intermediate- and high- 
resolution spectra of the quasar PKS 1448—232 
(zcm = 2.208; V = 16.9), observed at the VLT 
with the UVES spectrograph. The intermediate- 
resolution data have a spectral resolution of 
R w 45,000, corresponding to a velocity reso- 
lution of Av « 6.7 kms^i FWHM. These data 
were obtained and reduced as part of the ESO 
Large P rogramme "The Cosm ic Evolution of the 



IGM" (jBergeron et all 120021 ). The wavelength 



coverage of the intermediate resolution data is 
3050 - 10, 400 A. The S/N in the data varies be- 
tween 15 and 90 per spectral resolution element. 

The high- resolution data have R ^ 75, 000, 
corresponding to Aw « 4.0 kms~^ FWHM 
velocity resolution and were obtained with 
VLT/UVES in 2007 in an independent observing 
run (program ID 079.A-0303(A)). The wave- 
length coverage of the high resolution data is 
3000 - 6687 A. The raw data were reduced 
using the UVES pipeline implemented in the 
ESO-MIDAS software package. The pipeline re- 
duction includes flat-fielding, bias- and sky- 
subtraction and a relative wavelength calibra- 
tion. The individual spectra then have been cor- 
rected to vacuum wavelengths and coadded. The 
S/N in the high-resolution data is 20 — 70 per 
resolution element. 

2.2. Line-fitting method 

The detected absorption features that are asso- 
ciated with the two absorbers at Zabs = 2.1098 
and 2;abs — 2.1660 were fitted independently 
in both spectra (at intermediate resolution and 
high resolution) with Gaussian profiles using 
the CANDALF fitting routine 0, which uses a 
standard Levenberg-Marquard minimization al- 
gorithm. The program simultaneously fits the 
continuum and the absorption lines, delivering 
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Fig. 1: Absorption profiles for the Ovi absorber a,t z — 2.1098 in the high-resolution data (left 
panel) and the intermediate-resolution data (right panel). 



Table 1: Fit parameters for the system at z = 2.1098 







Z 




Ovi 


Civ 




Hi 






Ovi 


C IV 


Hi 


log[iV(cm-2)] 


b[kms^l] log[]V(cm"^)] 




log[JV{cm-2)] 


b [kms^l] 


lugh resolution data 


1 


2.10982 


2.10982 


2.10981 


14.27(±0.01) 


10.7(±0.2) 13.12(±0.01) 


7.5(±0.1) 


13.38(±0.04) 


19.6(±0.7) 


2 


2.11011 


2.11008 


2.11018 


13.50(±0.02) 


8.4(±0.4) 12.23(±0.03) 


6.1(±0.6) 


13.37(±0.04) 


28.6(±1.7) 


intermediate resolution data 


1 


2.10984 


2.10983 


2.10981 


14.32(±0.02) 


10.1(±0.2) 13.12(±0.01) 


7.1(±0.1) 


13.39(±0.01) 


20.6(±0.3) 


2 


2.11014 


2.11008 


2.11019 


13.49(±0.20) 


5.4(±0.5) 12.23(±0.03) 


5.3(±0.6) 


13.35(±0.01) 


26.4(±0.8) 



ion column densities, N, and Doppler parame- 
ters, 6, of for each absorption component. The 
continuum is modeled as a Legendre polynomial 
with an order of up to 4. The one-sigma fitting 



errors for N and b (as listed in Tables 1 — 3) 
are estimated using the diagonals of the Hesse 
matrix. 
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Fig. 2: Absorption profiles for the Ovi absorber at z = 2.1660 in the high-resolution data (left 
panel) and the intermediate- resolution data (right panel) . The strong absorption observed in O vi 
A1037.62 plot is a Sim line at z = 1.7236. 
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2.3. The Ovi system at z = 2.1098 

Fig. 1 shows the velocity profiles of O vi 
(AA1031.9, 1037.6), Civ (AA1548. 2, 1550.8), and 
Hi Lya and Ly/3 (AA1215. 7, 1025.7) for the 
z ~ 2.1098 absorber in the high-resolution data 
(left panel) and the intermediate-resolution data 
(right panel) . From the visual inspection of both 
panels we find no significant differences between 
the two data sets. The S/N ratios are better for 
the high resolution data, except for the C iv re- 
gion where this ratio is slightly better at interme- 
diate resolution. Therefore, the differences in the 
values for iV, 6, and z derived for the individual 
absorption components in the intermediate and 
high-resolution spectra are a result of the differ- 
ent S/N values in the two data sets. 

Two absorption components are detected in 
each of these ions. The O vi absorption is rela- 
tively strong compared to CiV. Hi absorption 
is weak compared t o other Ovi absorbe rs at 
similar redshift (e.g., Bergeron et al.|[2002l ) with 
a central absorption depth in the Hi Lya line 
of less than 70 percent. Note that the second, 
weaker, component of Hi Ly a and Ly^ absorp- 
tion associated with the high-ion absorption is 
blended, so that the true component structure 
of the Hi and the relative Hi column densi- 
ties and H i 6- values remain somewhat uncertain. 
The blending aspect is not taken into account in 
the formal error estimate for A'^ and h given in 
Table 1, which is based on the profile fitting. 
While for the stronger of these components the 
absorption of O VI, C IV, and H i is well aligned, 
there appears to be a small (< 10 kms~^) ve- 
locity shift between H i and the high ions in the 
weaker component (see Table 1). If real, this shift 
may indicate that the H i and the metal ions may 
not trace the same gas phase in the weaker ab- 
sorption component. Because of the blending of 
the Hi absorption, however, the reality of this 
shift remains unclear. 

We derive for the column densities, listed in 
Table 1, log 7V(0vi) « 14.3, log N{Gw) « 
13.1, and log A^(Hi) « 13.4 in the stronger 
of the two components. The resulting ion-to- 
hydrogen ratios of iV(0 Vi)/iV(Hi) ^ 8 and 



Table 2: Hi Ly a fit parameters for the absorp- 
tion system at z = 2.1660 





2 


log[Af(Hl) (cm^^)] 


b [kms-'] 


high resolution data 


1 


2.16466 


13.27 (±0.01) 


40.0 (±1.2) 


2 


2.16552 


14.49 (±0.07) 


29.9 (±1.2) 


3 


2.16605 


15.18 (±0.03) 


34.9 (±0.5) 


4 


2.16746 


13.16 (±0.03) 


19.2 (±0.6) 


5 


2.16765 


13.77 (±0.01) 


45.6 (±0.3) 




intermediate resolution 


data 


1 


2.16467 


13.21 (±0.02) 


37.4 (±1.3) 


2 


2.16559 


14.57 (±0.12) 


31.9 (±1.7) 


3 


2.16607 


15.13 (±0.05) 


35.5 (±0.8) 


4 


2.16749 


13.19 (±0.02) 


18.8 (±0.5) 


5 


2.16768 


13.76 (±0.01) 


46.8 (±0.3) 



7V(Civ)/iV(Hi) ~ 0.5 already indicate that the 
metallicity of this absorber mus t be fairly high 
(jBergeron fc Herbert-FortI l2005h . Note that be- 
cause of the blending problem in the Ly a and 
Ly /? lines the H i column density may be re- 
garded as upper limit, so that the ratios given 
above could be even higher. 



2.4. The Ovi system at z = 2.1660 

The Ovi system at z = 2.1660 exhibits a signif- 
icantly more complex absorption pattern than 
the absorber at 2 = 2.1098, as can be seen in the 
velocity profiles presented in Fig 2. Ovi absorp- 
tion is observed in eight individual absorption 
components, spanning a velocity range as large 
as ~ 300 kms~^. From the visual inspection it 
is further evident that the absorption pattern of 
O VI is different than those of the other detected 
intermediate and high ions (C III, C iv) and H i, 
although some of the components appear to be 
aligned in velocity space. As for the system at 
z — 2.1098, there are no significant differences in 
the absorption characteristics between the high- 
resolution data and the intermediate-resolution 
data. However, the S/N ratio is somewhat lower 
in the latter for the lines that are located in the 
blue part of the spectrum, so that the resulting 
fitting values for N, 6, and z for the individual 
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Table 3: Fit parameters for metal absorption in the absorber at ^ = 2.1660 







Z 




Ovi 




Civ 




Cm 






Ovi 


Civ 


Cm 


log[Ar(cm-2)] 


b [kms^l] 


log[Ar(cm-2)] 


blkms-'^] 


loglJV(cin-2)] 


6[kms-^] 


high resolution data 


1 


2.16518 






13.35(±0.03) 


15.4(±1.3) 










2 


2.16542 






13.02(±0.08) 


8.0(±1.2) 










3 


2.16669 


2.16569 


2.16566 


13.63(±0.02) 


13.6(±0.7) 


12.18(±0.0S) 


11.1(±1.7) 


12.68(±0.07) 


11. 1» 


4 


2.16600 


2.16600 


2.16602 


13.41(±0.03) 


9.9(±0.6) 


13.17(±0.05) 


9.3(±0.5) 


13.49(±0.05) 


15.1(±1.5) 


5 




2.16616 


2.16618 






12.96(±0.10) 


10.0(±1.6) 


12.78(±0.23) 


8.0(±3.7) 


6 


2.16638 


2.16638 


2.16633 


13.21(±0.04) 


13.7(±().8) 


12.53(±0.09) 


12.6(±2.2) 


12.22(±0.47) 


12.6" 


7 


2.16744 


2.16741 


2.16740 


12.75(±0.15) 


6.8(±0.8) 


11.90(±0.05) 


4.8'= 


12.17(±0.14) 


4.8* 


8 


2.16766 


2.16770 


2.16780 


13.01(±0.16) 


10.6(±0.4) 


11.84(±0.08) 


9.2'= 


12.32(±0.1B) 


9.2* 


9 


2.16796 


2.16792 


2.16789 


13.30(±0.08) 


10.9(±1.3) 


ll.B0(±0.12) 


4.8'= 


11.8B(±0.36) 


4.8=- 



intermediate resolution data 



1 


2.16521 






13.12{±0.03) 


8.3(±1.0) 










2 


2.16544 






13.13(±0.04) 


6.4(±0.9) 










3 


2.16600 


2.16567 


2.16557 


13.53(±0.05) 


10.6(±0.9) 


12.16(±0.05) 


11.1(±1.5) 


12.87{±0.11) 


14.4(±3.7) 
10.4" 


4 


2.16606 


2.16602 


2.16603 


13.73(±0.04) 


24.8(±3.1) 


13.23(±().()2) 


1().4(±().4) 


13.39(±0.08) 


6 




2.16619 


2.16626 






12.73{±().()!)) 


7.5(±0.9) 


13.12(±0.15) 


7.5" 


6 


2.16646 


2.16637 


2.16635 


12.97(±().ll) 


10.4(±2.0) 


12.(i7(±().()8) 


15.3(±2.6) 






7 


2.16748 


2.16742 




12.98(±0.03) 


7.3(±0.9) 


11.67(±0.17) 


4.8(±2.3) 






8 


2.16767 


2.16758 




12.93(±0.05) 


2.8(±0.9) 


12.24(±0.07) 


22.2(±3.1) 






9 


2.16796 


2.16797 




13.50(±0.01) 


12.7(±0.6) 


11.85(±0.06) 


4.8(±1.2) 







Fixed to 6c III — ^CIV 



Fixed to b-valuc derived from the intermediate resolution data 
*^ Lower limit, fixed to the minimal value 



absorption components differ slightly (Tables 2 
and 3). 

We have modeled the Hi absorption by si- 
multaneously fitting Lya and Ly/3 in four ab- 
sorption components (components 2 — 5; see 
Table 2), obtaining column densities between 
13.2 < log iV(Hi)< 15.2. One additional com- 
ponent (component 1) is present in the Lya ab- 
sorption, but is blended in Ly /3 (sec Fig. 2), so 
that A^(Hi) was derived solely from Lya. Note 
that for the H i fit we have not tried to tie the H i 
component structure to the structure seen in the 
the metal ions, as this requires knowledge about 
the physical conditions in the absorber. This as- 
pect will be discussed in detail in Sect. 4.2, where 
we try to reconstruct the H i absorption pattern 
based on a photoionization model. Instead, we 
have fitted the H i absorption with the minimum 
number of absorption components required to 
match the observations (Fig. 2, lowest panel) and 
to obtain an estimate on the total H I column in 
the absorber. 

By summing over the column densities in 
the individual absorption components we de- 
rive total column densities of log A^(0 vi) « 
14.2, log 7V(Ciii) « 13.7, log 7V(Civ) « 13.5, 



and log iV(IIi) w 15.3. The resulting ion-to- 
hydrogen ratios of iV(0 vi)/7V(Hi) ~ 0.1 and 
N{Civ)/N{m) ~ 0.02 (representing the av- 
erage over all components) are substantially 
smaller than in the z — 2.1098 system, point- 
ing toward a lower (mean) metallicity of the ab- 
sorber. 

The complexity in the absorption pattern of 
the various species in this system and the large 
velocity spread suggests that this absorber arises 
in an extended multi-phase gas structure. 

3. Ionization modeling and physical 
conditions in the gas 

To infer information on the physical properties of 
the two Ovi absorbers towards PKS 1448-232 
we have modeled in detail the ionization condi- 
tions in these systems. Since the two absorbers 
at z = 2.1098 and z = 2.1660 have redshifts 
close to the quasar redshift (zqso = 2.208), it is 
necessary to check whether the two systems lie 
in the proximity zone of the background quasar 
and are influenced by its ionizing radiation. 

With the above given redshifts, the two 
absorbers have velocity separations from the 
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Table 4: Modelled column densities for the absorber at z = 2.1098 







log 


[N (cm" 
















V [kms ^\ 


Civ 


Ovi 


Hi 


log [tih (cm"^)] 


log Z 


log [r(K)] 


L [kpc] 


./hi 


1 





13.12 


14.27 


13.38 


-4.20 


-0.24 


4.54 


19.9 


-5.21 


2 


+25 


12.23 


13.50 


13.37 


-4.25 


-1.02 


4.64 


30.5 


-5.35 


2 


+25 


12.23 


13.50 


12. 57=^ 


-4.28 


-0.24 


4.57 


4.7 


-5.32 



^ Our best H I guess in the model for the second component with fixed metallicity 



QSO of (5w2.io98 ~ 9000 kms^^ and 5v2.ima ~ 
4000 kms^^ and thus the absorber at z = 
2.1660 can be regarded (depending on the 
definition) as an associated system. With a 
(monochromatic) luminosity at the Lyman Limit 
of ^912 = 3.39 X lO^i ergs-iHz-\ the size of 
the sphere-of-influence of the ionizing radiation 
from PKS 1448-232 is known to be 6.7 Mpc, cor- 
respond ing to a velocit y separation of ~ 1400 
kms-i (|Fox et al.ll2008[) . Therefore, it is safe to 
assume that the ionizing radiation coming from 
PKS 1448—232 itself has no measurable influence 
on the ionization conditions in the two O vi sys- 
tems. 

The small &- values measured for O VI, C iv, 
and H i indicate that coUisional ionization is not 
the origin for the presence of O vi in the gas. It 
is common to assume that the observed Doppler 
parameters (6obs) are composed of a thermal and 
a turbulent component (6th and &turb, respec- 
tively), so that 6oi,g = + 6turb- The thermal 
component can be expressed by — 2kT/m, 
where T is the gas temperature and m is the 
mass of the considered ion. The Doppler param- 
eters measured for the O VI components in the 
two absorbers are all 6 < 16 kms~^ and many 
of them are 6 < 10 kms~^ (see Tables 1 and 
3), indicating that T < 10^ K. This value is 
below the peak temperature of Ovi in a colli- 
sional ionization equilibriu m (T ~ 3 x 10^ K; 
Sutherland fc Dopital 1993 ) ; it is also lower than 
the temperature range expected for O vi arising 
in turbulent mixing layers in the interface re- 
gio ns between cold and ho t gas (T = 10^ — 10^ 
K; iKwak fc SheltonI |2010[ ). Consequently, pho- 
toionization by the hard UV background remains 
as the only plausible origin for the presence 



of O VI in the two high-ion absorbers towards 
PKS 1448-232. 

Based on these considerations, we have mod- 
eled the ion column densities in the two O vi sys- 
tems using th e photoionization c ode CLOUDY 



(version C08: iFerland et al.llT998l ). For this 



we 

have assumed a solar relative abundance pattern 
of O and C and an optically thin plane-parallel 
geometry i n photoionization equilibrium, ex- 
posed to a iHaardt fc Madaul (200 Ij UV back- 
ground spectrum at z = 2. 16, which is norm al- 
ized to log J912 = -21.15 (|Scott et al.ll2000[ ) at 
the Lyman limit. 

We further assume that each of the observed 
velocity components is produced by a "cloud", 
which we model as an individual entity. As in- 
put parameters we consider the measured col- 
umn densities of Cm (only for the z = 2.1660 
absorber). Civ, Ovi, the metallicity Z (in so- 
lar units), and the hydrogen particle density Uh- 
The metallicity of each cloud and the hydrogen 
density were varied in a range appropriate for 
intergalactic clouds (i.e., —3 < log Z < and 
-5 <log riH < 0). 

We then applied the following iterative mod- 
eling procedure. In a first step, CLOUDY was 
run with a set of values of Z, uh and A^(Hi), 
where iV(Hl) is constrained by the observations. 
In a second step, the corresponding values of 
iV(Ciii), 7V(Civ), and iV(Ovi) were calculated. 
The output was compared with the observed col- 
umn densities and, in case of a mismatch, the in- 
put parameters Z and uh were adjusted for the 
next iteration step. This process was repeated 
until the differences between the output column 
densities and the observed values became negligi- 
ble and we obtained a unique solution. Next to 
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the ion column densities, our CLOUDY model 

provides information on the neutral hydrogen 
fraction, /hi, the gas temperature, T, and the 
absorption pathlength, L = N{'H.i)/{fm'nu)- 

3.1. The system at z = 2.1098 

As mention earlier, absorption by O VI and C iv 

is well aligned in both components in this sys- 
tem, while the true component structure of the 
H I is uncertain because of blending effects in the 
Ly a and Ly /3 lines. Because of the alignment of 
Ovi and Civ we assume a single-phase model, 
in which each of the two components (clouds) 
at u = and +25 kms'^ in the z = 2.1098 
rest frame hosts Ovi, CiV, and Hi at column 
densities similar to the ones derived from the 
profile fitting. Consequently, we have chosen log 
7V(Hi) = 13.37 and 13.38 as input for the 
CLOUDY modeling and followed the procedure 
outlined above. The results from the CLOUDY 
modeling of the z = 2.1098 absorber are summa- 
rized in Table 4. Our model reproduces well the 
observed O vi and C iv cohimn densities in both 
components, if the clouds have a density of log 
riH ~ ^4.2, a temperatmc of log T w 4.6, and 
a neutral hydrogen fraction of log /hi ~ —5.3. 
However, to match the observations, the second 
component (at +25 kms~^) in our initial model 
(Table 4, first two rows) needs to have metallic- 
ity of log Z — —1.02, which is ^ 0.8 dcx lower 
than for the other component (log Z = —0.24). 
The absorption path lengths are ^ 20 kpc for 
the component at kms~^ and ~ 30 kpc for the 
component at +25 kms~^. 

Because of the blending problem in the Hi 
Lya and Ly/3 absorption, which affects par- 
ticularly the estimate for A'' (Hi) in the cloud 
at +25 kms"-'^ (Fig. 1), we have set up a sec- 
ond CLOUDY model in which we have tied 
the metallicity of the +25 kms~^ component 
to the metallicity of the other component (log 
Z = -0.24), but leaving iV(Hi) for this compo- 
nent as a free parameter. ^Prom this we derive a 
value of log Af(Hi) = 12.57 for the cloud at +25 
kms~^ and the absorption path-length reduces 



to L = 4.7 kpc. In view of the blending, we re- 
gard this model as more plausible compared to 
the model with two different metallicities and a 
larger absorption path-length. 

Summarizing, our CLOUDY modeling sug- 
gests that the z = 2.1098 absorber towards 
PKS 1448— 232 represents at relatively simple, 
metal-rich O vi absorber in which the high ions 
O VI and C iv coexist in a single gas-phase. 

3.2. The system at z = 2.1660 

We started to model this system with CLOUDY, 
again under the assumption of a single gas-phase 

hosting the observed intermediate and high ions 
Cm, Civ, and Ovi in the various subcompo- 
nents. However, during the modeling process it 
quickly turned out that it is impossible to match 
the observed column densities of C ill and O VI 
in a single gas-phase in the components, where 
these two ions are aligned in velocity space. Our 
modeling instead indicates that the C ill absorp- 
tion must arise in an environment that has a 
relatively high gas density and that is spatially 
distinct from the O VI phase. In a second step, we 
have tried to tie together the high ions C IV and 
Ovi in a single gas phase (as for the z = 2.1098 
system) in the relevant absorption components, 
ignoring the Cm phase. Again, this approach 
does not deliver satisfying results, as we obtain 
for some of the components, for which C iv/0 VI 
is constrained by observations, very low gas den- 
sities and very large absorption patlilengths on 
Mpc scales, which are highly unrealistic. Given 
the fact, that the overall component structure of 
Ovi and Civ are substantially different in this 
system (Fig. 2), this result is not really surpris- 
ing. 

The only modeling approach, for which we 
obtain realistic results on gas densities, tempera- 
tures and absorption path-lengths in this system 
and its subcomponents is a two-phase model, in 
which C III coexists with C IV and part of the H I 
in one (spatially relatively confined) phase, and 
O VI and part of the H I in a second (spatially rel- 
atively extended) phase. The coexistence of C ill 
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Table 5: Modelled column densities for the Ciii/Civ absorbing phase in the z = 2.1660 absorber 





V [km s ^ ] 


Cm 


log {N ( 
Civ 


cm-^)] 
O VI 


Hi" 


log [riH (cm ^)] 


log Z 


log [T(K)] 


L [kpc] 


/h 


3 


-28 


12.68 


12.18 


10.97 


14.51 


-2.74 


-1.7 


4.42 


0.3 


-3.68 


4 


+0 


13.49 


13.17 


12.25 


14.18 


-2.97 


-1.7 


4.46 


4.1 


-3.95 


5 


+ 16 


12.78 


12.96 




14.51 


-3.56 


-1.7 


4.58 


16.3 


-4.63 


6 


+37 


12.22 


12.53 


12.86 


14.08 


-3.71 


-1.7 


4.61 


12.3 


-4.79 


7 


+ 134 


12.17 


11.90 


10.92 


13.26 


-2.93 


-1.0 


4.38 


0.04 


-3.84 


8 


+ 162 


12.32 


11.84 


10.49 


13.57 


-2.67 


-1.0 


4.34 


0.02 


-3.55 


9 


+ 182 


11.85 


11.50 


10.38 


12.99 


-2.84 


-1.0 


4.37 


0.01 


-3.73 


* Our 


best H I guess 


in the 


models. 

















and C IV in one phase is further suggested by the 
fact that C ill and C iv absorption is well aligned 
in velocity space (see Fig. 2). The results from 
this two-phase model are presented in Tables 5 
and 6. A critical issue for the modeling of this 
complex multi-phase absorber with its many ab- 
sorption components is an assumption for the 
neutral gas column density in each subcompo- 
nent (and phase). Since in the Hi Lya and Ly/3 
absorption most subcomponents are smeared to- 
gether to one large absorption trough, the obser- 
vational data provide little information on the 
distribution of the Hi column densities among 
the individual components. Yet, the data give a 
solid estimate for the total H i column density in 
the absorber (log N k. 15.3; see Sect. 2.4), which 
must match the sum of A^(Hi) over all subcom- 
ponents considered in our model. Consequently, 
we included in our iteration procedure the con- 
straints on A^(Hi)tot and the shape of the (to- 
tal) H I absorption profile. The latter aspect also 
concerns the choice of the gas temperature in 
the model, as T regulates the thermal Doppler- 
broadening and thus the width of the modeled 
Hi lines. We have modeled the Hi width follow- 



ing the aproach of iDing et al.l (|2003| ). 



With these various constraints we first mod- 
eled the Ciii/Civ phase in the absorber. 
However, because of the extremely complex pa- 
rameter space, we did not find a unique solution 
for (T, riH, among the individual components, 
but had to make further constraints. Since the 
individual components observed in C iii/C iv are 
very close together in velocity space, we assumed 
they all have the same metallicity and, based on 



the Z range allowed in the model, we set log 
Z = —1.5 for all subcomponents. This model 
was able to match the observed column densi- 
ties of these two ions in the individual subcom- 
ponents, but did not match well the gross shape 
of the overall Hi absorption, suggesting that 
the metallicity in this absorber is non-uniform 
among the individual absorption components. 
Therefore, we refined our model, now using two 
different metallicities: log Z = —1.7 for the sat- 
urated Hi components and log Z = —1.0 for 
the weaker H i components (see Tables 5 and 6 
for details). Although not perfect, this model de- 
livers a satisfying match between the modeled 
spectrum and the UVES data. 



Adopting this model, we find that the 
Ciii/Civ absorbing components have tempera- 
tures between log T = 4.3 and 4.6, densities be- 
tween log nn = —3.7 and —2.7, and neutral gas 
fractions between log /hi = —4.8 and —3.6 (see 
Table 5). The absorption path lengths vary be- 
tween 0.3 and 16.3 kpc for the components with 
log Z = —1.7, and between 0.01 and 0.04 kpc 
for the components with log Z = —1.0. These 
numbers suggest that the Ciii/Civ absorbing 
phase resides in relatively small and confined gas 
clumps. This scenario is supported by the small 
turbulent 6- values of < 6 kms^^ for the subcom- 
ponents that we derive in our model. Note that 
in Table 5 we also list the predicted column den- 
sities for O VI, which are typically 1 — 2 orders 
of magnitude below the observed ones in this 
absorber. This, again, underlines that Ciii/Civ 
and O VI must reside in different gas phases with 



N. Draganova, P. Richter, and C. Fechner: O vi systems towards PKS 1448—232 



11 



different physical conditions to explain the ob- 
served column densities. 



Finally, we have modeled the O vi absorbing 
phase in the z — 2.1660 absorber, based on the 
observed Ovi column densities. Since there are 
no ions other than H i and O vi that could pro- 
vide information about the physical conditions 
in this phase, we fixed the metallicity of the gas 
to log Z = —1.7 and log Z = —1.0 (equal to the 
Ciii/Civ phase) and constrained the tempera- 
ture range [Tmin, Tmax] in the CLOUDY models 
based on the observed line widths of O vi (giving 
Tmax) and the modeling results of the Ciii/Civ 
phase (giving Tmin for all components except the 
first two). The results of this model are shown in 
Table 6. We derive gas densities in the range log 
riH = —4.6 to —3.2 and neutral gas fractions in 
the range log /hi = —5.8 to —4.6. The absorp- 
tion path length varies between 19.8 and 83.3 
kpc for the components with log Z — —1.7, and 
between 1.3 and 38.3 kpc for the ones with log 
Z = —1.0. The mismatch in N{On\) between 
the model and the data for components one and 
nine (see Table 6) points towards a metallic- 
ity distribution among the individual absorption 
components that is even more complex than the 
one assumed in our model. Despite this (mi- 
nor) concern, our CLOUDY modeling for O vi 
provides clear evidence that the O vi absorbing 
phase has substantially lower gas densities than 
the C iii/C IV absorbing phase and is spatially 
more extended. 



In summary, our CLOUDY modeling of the 
z — 2.1660 absorber suggests that this system 
represents a complex multi-phase gas structure, 
in which a number of cooler, C iii/C iv absorbing 
cloudlets are embedded in a spatially more ex- 
tended, O VI absorbing gas phase spanning a to- 
tal velocity range of ~ 300 kms""'^. Although the 
metallicity is not well constrained in our model, 
it suggests that log Z < —1 in the absorber, 
which is ^ 0.8 dex below the value obtained for 
the system at z — 2.1098. 



4. Discussion 

Our detailed analysis of the two O vi absorbers 
at z = 2.1098 and z = 2.1660 towards the quasar 
PKS 1448-232 displays the large diversity and 
complexity of high-ion absorbers at high red- 
shift. 

During the past years, a number of 
studies using both optical observations 



(e.g., Bereeron et al. 


2002; Carsweh et al. 


2002; Simcoe et all 


2002L 20041 20061 


Bergeron & Herbert-Fort 


2005; Affuirre et al. 


2008) and numerical simulations (e.g., 


Faneanoetal. 2007; Kawata & RauchI 2007) 


have been dedicated to 


investigate the prop- 



erties of high-redshift Ovi systems and their 
relation to galaxies. 

Based on their survey of Ovi absorbers 

in the redshift range z = 2.0 — 2.6, 

Bergeron fc Herbert-Fori] (|2005h suggested that 



Ovi systems may be classified into two dif- 
ferent populations: metal-rich absorbers ("type 
1") that have large A^(0 Vi)/A^(H i) ratios and 
that appear to be linked to galaxies and galac- 
tic winds, and metal-poor absorbers ("type 0") 
with small A^(0 vi)/A^(Hi) ratios that trace 
the intergalactic medium. The two absorbers 
towards PKS 1448—232 discussed in this pa- 
per do not match the classifica tion scheme of 
Bergeron fc Herbert-Fod] (|2005l ). The absorber 



at z = 2.1098 has a very large iV(0 vi)/Af(H i) 
ratio of ~ 8 (i.e., it is of type 1); it is a simple, 
single-phase, metal-rich system with a metal- 
licity slightly below the solar value. Yet, this 
system is completely isolated with no strong 
Hi Lya absorption within 1000 kms^^. In 
contrast, the absorber at z = 2.1660 has a 
N{Ovi)/N(Hl) ratio of only ~ 0.1 and a 
metallicity of 0.1 sol ar or lower (i.e., it is of 



type according to iBergeron fc: Herbert-Fort 
20051 ) ■ However, this absorber is a complex 
multi-phase system with a non-uniform metal- 
licity, suggesting that originates in a circum- 
gala ctic environment. While this mis match with 
the iBergeron" fc Herbert-PbT^ (|2005l) classifica- 
tion scheme certainly has no statistical relevance 
for the general interpretation of Ovi absorbers 
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Table 6: Modelled column densities for the Ovi absorbing phase in the z = 2.1660 absorber 





V [km s ^\ 


log [N ( 
Ovi 


cm '^)\ 
Hi" 


log [riH 


(cm-3)] 


log Z 


log [T(K)] 


L [kpc] 




/h 


1 


-78 


12. 96'^ 


13.17 


< - 


3.76 


-1.7 


<5.36 




<19.8 


> 


-5.86 


2 


-55 


13.02 


14.10 


< - 


3.77 


-1.7 


<4.79 




<25.8 


> 


-5.03 


3 


-30 


13.63 


14.51 


-3.91 


.. -3.31 


-1.7 


4.42 ... 5.24 


58.4 . 


.. 67.7 


-5.50 


... -4.84 


4 


+0 


13.41 


14.18 


-3.99 


.. -3.88 


-1.7 


4.46 ... 4.98 


42.6 . 


.. 83.3 


-5.36 


... -4.95 


6 


+36 


13.21 


14.08 


-3.92 


.. -3.23 


-1.7 


4.61 ... 5.26 


32.1 . 


.. 21.5 


-5.52 


... -5.00 


7 


+ 136 


12.75 


13.26 


-3.72 


.. -3.72 


-1.0 


4.38 ... 4.64 


1.3 


... 2.0 


-4.82 


... -4.61 


8 


+ 157 


13.01 


13.57 


-3.70 


.. -3.52 


-1.0 


4.34 ... 5.03 


2.2 


... 6.1 


-5.19 


... -4.56 


9 


+ 185 


13.25= 


12.99 


-4.60 


.. -4.25 


-1.0 


4.37 ... 5.06 


32.3 . 


.. 38.3 


-5.76 


... -5.48 



" Our best H I guess in the models 

Observed log JV(Ovi) = 13.35 
" Observed log JV(0 Vl) = 13.30 



at high redshift, the results suggest that for a 
thorough understanding of highly-ionized gas at 
high redshift the absorption characteristics of 
Ovi systems may be too diverse for a simple 
classification scheme based solely on observed 
(and partly averaged) column density ratios of 
O VI, Hi and other ions. 

One critical drawback of many previous O vi 
surveys at high z is that they often consider only 
simplified models for the ionization conditions 
in their sample of high-ion absorbers, so that 
the multi-phase character of the gas and pos- 
sible ionization conditions far from a photoion- 
ization equilibrium are only insufficiently taken 
into account. As pointed out by Fox (201 1), 
single-phase, single-component ionization mod- 
els, if applied, will deliver physically irrelevant 
results for most of the Ovi systems at high 
z. This implies that previous estimates of the 
baryon- and metal-content of O VI absorbers at 
low and high z possibly are afflicted with large 
systematic uncertainties. 

One firm conclusion from many previous ob- 
servational and theoretical studies of high-ions 
absorbers is that a considerable fraction of the 
O VI systems at low and high z must arise in the 
metal-enriched circumgalacti c environment of 
(star -f orming) galaxies (e.g., Wakker fc Savagel 



2009: Prochaska et al. 2011; Fox et al 



Tepper-Garcia et al..,2011; jFangano et al 



20111: 



20071) 



Thus, the complex absorption pattern observed 
in the z = 2.1660 system towards PKS 1448-232 
and many other O vi absorbers at high z may 



reflect the complex gas distribution of enriched 
gaseous material that was ejected from galax- 
ies into th e IGM during their wind-blowing 



phas e (e.g., Kawata, h Rauchll2007l ). In this con- 



text, Schave et al. ( 2007f ) suggested that the in- 



tergalactic metals have been transported from 
galaxies through galactic winds and reside in 
the form of dense, small and high-metallicity 
patches within large hydrogen clouds. These au- 
thors point out that much of the scatter in the 
metallicities derived for high-redshift absorbers 
could be explained by the spatially varying num- 
ber of the metal-rich patches and the different 
absorption path lengths through the surround- 
ing metal-poor intergalactic filament instead of 
an overall (large-scale) metallicity scatter in the 
IGM. In this scenario, the substantial differ- 
ences in the metallicities of the two O vi systems 
towards PKS 1448—232, and even the intrinsic 
metallicity variations within the z = 2.1660 sys- 
tem, could be explained by the different geome- 
tries of the absorbing structures, suggesting that 
much of the H i that is associated with the metal 
absorption in velocity space, arises in a spatially 
dis tinct region. A similar conc lusion was drawn 
bv iTepper-Garda et al.l ( 2011 ). who studied the 
nature of Ovi absorbers at low z using a set 
of cosmological simulations. Note that also ab- 
sorbers with larger Hi column densities, such as 
Lyman-limit systems (LLS) and damped-Lyman 
a systems (DLAs) , sometimes exhibit abundance 
variations am ong the different velocity subcom- 
ponents (e.g.. lRichter et all2005 : Prochter et al 
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This indicates that the metals in the gas 
surrounding high z galaxies are not well mixed. 

The observed velocity differences between 
O VI and other ions and the multi-phase nature 
of the gas provide further evidence for an inho- 
mogeneous metallicity and density distribution 
in intervening high-ion absorbers. It is an inter- 
esting fact that the velocity misalignment ap- 
pears to concern only the O VI absorbing phase 
in high-ion absorbers at high redshift, while 
other high ions such as N v and C IV generally 
appear to be well aligned with H i, even in sys- 
tems that exhibit a complex velocity -component 
structure (jPechner fc Richten 120091 ). This puz- 
zling aspect underlines that additional detailed 
studies of individual O VI absorption systems 
could be very important for our understanding 
of intergalactic and circumgalactic gas at high 
redshift, as this ion traces a metal-enriched gas 
phase that cannot be observed by other means. 

5. Summary and outlook 

In this paper, we have investigated two O VI 
absorbers at z = 2.1098 and z = 2.1660 to- 
wards the quasar PKS 1448-232. For this, we 
have used high- {R ^ 75, 000) and intermediate- 
resolution {R w 45, 000) optical spectra obtained 
with the VLT/UVES instrument and CLOUDY 
photoionization models. 

The Ovi system at 2: = 2.1098 is character- 
ized by strong O VI absorption and weak H i ab- 
sorption in a relatively simple, two-component 
absorption pattern. The absorption by O vi, 
C IV, and H i are well aligned in velocity space, 
indicating that they trace the same gas phase. 
From a detailed photoionization modeling of this 
system we derive a metallicity of 0.6 solar, a 
characteristic density of log rin ~ —4.2, a tem- 
perature of log T fa 4.6, and a total absorption 
path length of ^ 30 kpc. The absorber is isolated 
with no strong H I Ly a absorption within 1000 
kms~^. 

The Ovi absorber at z = 2.1660 repre- 
sents a complicated, multi-component absorp- 
tion system with eight relatively weak and nar- 



row O VI absorption components spanning al- 
most 300 kms~^ in radial velocity. The Ovi 
components are accompanied by strong Hi ab- 
sorption and C III, C IV absorption. The O vi 
component structure differs from that of H i 
and C IV, indicating a multi-phase nature of the 
absorber. Our photoionization modeling with 
CLOUDY suggests the presence of (at least) two 
distinct gas phases in this system. Cm, Civ 
and most of the Hi appear to coexist in sev- 
eral, relative compact cloudlets at gas densities 
of log riH « —3.7 to —2.7, temperatures of log 
T w 4.3 — 4.6 and absorption path lengths of 
< 16 kpc. O VI appears to reside in a highly ion- 
ized, more extended gas phase at densities in the 
range log nn ~ —4.6 to —3.2, temperatures be- 
tween log T Ki 4.3 and 5.3, and absorption path 
lengths up to 83 kpc. While the exact metal- 
licity of the absorber is not well constrained, 
our modeling favours a non-uniform metal abun- 
dance among the individual absorption compo- 
nents with (at least) two different metallicities 
of log Z = -1.7 and log Z = -1.0. 

Our study displays the large diversity and 
complexity of Ovi systems at high redshift. 
We speculate that some of the observed differ- 
ences between the two high-ion absorbers to- 
wards PKS 1448—232 could be a result of a inho- 
mogeneous metallicity and density distribution 
in the photoionized IGM. Our study indicates 
that multi-phase, multi-component high-ion ab- 
sorbers like the one at 2; = 2.1660 demand a 
detailed ionization modeling of the various sub- 
components to obtain reliable information on 
physical conditions and metal-abundances in the 
gas. We conclude that a rather large effort is re- 
quired to achieve a more complete view on the 
nature of O vi absorbers at high redshift. 

For the future, we are planning to continue 
our investigation on these systems by using a 
larger sample of O vi absorbers in high-quality 
UVES archival data and compare their absorp- 
tion characteristics with artificial spectra gener- 
ated from numerical simulations of star-forming 
galaxies and their intergalactic environment. 
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